double torad(double deg) {return deg / 180 * pi;} 
Point3 get_coord(double R, double lat, double lng) {
    lat = torad(lat); lng = torad(lng);
    return P(R * cos(lat) * cos(lng), R * cos(lat) * sin(lng), R * sin(lat));
}
